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ABSTRACT 

We show that the collapsar model of gamma-ray bursts results in a series of 
successive shocks and rarefaction waves propagating in the "cork" of stellar material 
being pushed ahead of the jet, as it emerges from the massive stellar progenitor. 
Our results are derived from analytical calculations assuming a hot, ultrarelativistic 
1-D flow with an initial Lorentz factor Tj ~ 100. The shocks result in a series of 
characteristic, increasingly shorter and harder thermal X-ray pulses, as well as a 
non-thermal 7-ray pulse, which precede the usual non-thermal MeV 7-rays. We 
consider jets escaping from both compact (CO or He dwarf) and blue supergiant stellar 
progenitors. The period, fluence and hardness of the pulses serves as a diagnostic for 
the size and density of the outer envelope of the progenitor star. 



Subject headings: Gamma-rays: bursts - X-rays: bursts - Cosmology: miscellaneous 
Stars: supernovae: general 



1. Introduction 



The non-thermal radiation from shocks in a relativistic fireball jet after emerging from a 
collapsing massive stellar system (collapsar) is the leading theoretical explanation for the family of 
"long" gamma-ray bursts (GRB) with 7-ray durations longer than about i& ~ 10 s (Woosley 1993, 
Paczyhski 1998; see Meszaros 2002 for a review). The jet is thought to arise from a brief accretion 
episode of a rotating debris torus around a central black hole formed as the stellar core collapses. 
The nature and history of the massive star, however, is so far largely a matter of hypothesis 
and calculation, and alternative stellar progenitor systems are conceivable, both for the long and 
especially for the shorter bursts. 

The duration of the TeV neutrino burst expected from shocks in the jet before it emerges is 
a possible diagnostic for the pre-emergence jet history and for the dimensions or column density 
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of the progenitor stellar system (Meszaros & Waxman 2001). Previous analytical work on shock 
and jet emergence has been done by Matzner & McKee 1999, Meszaros & Rees 2001, Matzner 
2002, and numerical calculations of relativistic jets have been done by Marti, Miiller, Font, Ibahez 
& Marquina, 1997; Aloy, Miiller, Ibahez, Marti & MacFadyen, 2000; and Zhang, Woosley & 
McFadyen, 2002. Here we concentrate on possible photon signatures of the jet as it emerges, which 
precedes the usual gamma-ray emission and the subsequent longer wavelength afterglow. This 
has been considered by Ramirez- Ruiz, McFadyen & Lazzati (2001), who infer an X-ray precursor 
to the burst from the portion of the stellar envelope shocked by the jet. Here we investigate in 
some detail, using analytical 1-D methods in the ultrarelativistic limit, the jet emergence process 
and the shock heating and expansion of the plug of stellar material propelled by it, resembling an 
ejected cork. 

In the initial stage of ejection the cork is optically very thick to scattering, and it experiences 
successive shocks and rarefactions as the optical depth decreases during its overall expansion. This 
leads to an X-ray photon precursor to the usual GRB emission, consisting of a pattern with an 
initial rise and decay followed by a brighter and harder peak, possibly repeated a few times with 
decreasing period. The characteristic timescale of the initial decay and the ratio of its amplitude 
to that of the peak differ substantially depending on the density and extent of the envelope, thus 
providing a potential diagnostic for the pre-burst stellar configuration. These precursor patterns 
should be detectable with instruments to be flown in the next few years. 



2. Jet propagation 

We assume a jet with isotropic equivalent luminosity Li so <; 10 52 Ls2 erg s _1 and Lorentz 
factor Tj ~ 10 2 . As the jet advances through the star, it drives forward a bow shock into the star. 
The jet is capped by a "cork" composed of a termination shock advancing into the star, and a 
reverse shock moving back into the jet, where the tenuous and highly relativistic jet gas is suddenly 
decelerated and compressed. Both the shocked jet plasma and the shocked stellar plasma advance 
inside the star with a jet head Lorentz factor <C Tj, corresponding to a velocity (3h. Figure [j] 
presents a schematic description of the multiple shock structure. In the rest frame of the shocked, 
decelerated plasma, the fast jet moves with a Lorentz factor ~ Tj/2Th. During propagation in the 
He core, the high density ahead of the jet, p ~ 1 g cm -3 , implies deceleration of the jet plasma 
to a sub-relativistic velocity. However, beyond the C/O or He core edge at r Z 10 10 - 10 11 cm, 
there is a sharp drop in the density, e.g. the density of an H envelope, if present, is typically 
p ~ 10 -5 — 10 _7 g cm -3 , and the jet can accelerate to relativistic velocity. We have used Matzner 
& McKee's parametrization of the density drop near the edge of the star, r = R*, 

where n = 3(3/2) for radiative (convective) envelope. Eq. (||) holds for a polytropic plasma at 
radii beyond which the envelope contribution to the gravitational potential is small. We consider 
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three progenitor examples: i) the edge of the C/O progenitor model of SN1998bw (Woosley, 
Eastman & Schmidt 1999), approximately described by = 5 x 10 2 g/cm 3 , ii* = 10 10 cm and 
n = 3; ii) the outer part of the He mantle of the pre-supernova model, evolved by Woosley, 
Langer & Weaver (1993) from a 35M Q main sequence star, approximately described by Eq. (|l]) 
with = 2g/cm 3 , i?* = 10 11 cm and n = 3; and hi) the envelope of the blue supergiant (BSG) 
progenitor models of SN1987A (Shigeyama & Nomoto 1990,Arnett 1991), approximately described 
by p* = 3 x 10 _5 g/cm 3 , R* = 3 x 10 12 cm and n = 3. Taking in Eq.(|l|) a radius r ~ 0.9i?* we take 
as approximate reference quantities the values p ~ [1, 10~ 3 , 10~ 7 ]g cm -3 at r ~ [10 10 , 10 11 , 10 12 ' 5 ] 
cm for a C/O core, He core and BSG star H envelope respectively. 

The pressure behind the forward shock is given, for rj, > 1, by = (4/3)r 2 l pc 2 , and for 
flh <C 1, by Ph = [(7 + l)/2]p/3^c 2 , where 7 = 4/3 is the adiabatic index of the plasma (radiation 
dominated). Equating Ph with the pressure behind the reverse shock, P r = (A/3)(Tj/2Th) 2 njm p c 2 
for Th S> 1 and P r = (4/3)T 2 njm p c 2 for Ph "C 1, where the jet proper proton density 
rij = Li so / (4Trr 2 T 2 m p c 3 ) (For a discussion of relativistic flows and shock jump conditions, see e.g., 
Begelman, Blandford & Rees 1984). The extreme relativistic limit gives (Meszaros & Waxman 
2001) 

r 1 fi ri/4/ 1/2 1/4 tt. 

fc 1 nor 1 ^ / 1/2 1/4 w V^J 
I 0-9L 5 ^ / ri [ pl % He 

for the H-envelope of a BSG (H) or a He star, while the non-relativistic limit gives 

~ / L6L 52 2 /ni/5y 3 2 He; r 
W i 1/2, 1/2 

1 0.6L^ 2 /noPo /2 CO 

for the He or CO cores, where more exactly (3h < 1 and > 1 would hold; Lj so = 10 52 L52erg s _1 
is the jet isotropic equivalent luminosity, r = W x r x is the shock radius, and p = W x p x g cm -3 is 
the corresponding density at that radius. These values are for the main part of the cork, including 
the shocked jet and the bulk of the shocked stellar material. Thus, for r ~ 0.9-R* in all cases 
the jet head velocity is at least marginally relativistic, and this approximation gets better as 
r — > i?*. This implies that, even though the jet injected at the base of the jet is assumed to be 
ultrarelativistic, Tj = 10 2 T2 3> 1, the jet head advances sub-relativistically for most of its initial 
crossing of the stellar interior. This is due to the initial inertia of the dense stellar material which 
needs to be pushed ahead and aside (e.g. Meszaros &; Rees, 200 for an analytical discussion). A 
similar relativistic injection and subrelativistic jet head advance occurs in AGN jets encountering 
the intergalactic medium, e.g. Begelman et aZ.1984. The jet crossing time out to the He (or 
CO) core edge is of order t core ^ 30rn s for a jet which has reached Vh ~ c near the core edge, 
suggesting that compact (r <, 10 11 cm) stars such as He stars or pre-WR stars are likely GRB 
progenitors (see also MacFadyen et al. 2001 and Matzner 2002). 

However, if there is an H-envelope beyond the C/O/He core, where the density drops steeply 
by 10 -5 — 10 -7 , momentum conservation indicates that Yh for the entire cork (equation [^]) can 
become significantly larger at radii just beyond r core = 10 n i?n cm. This jump in Th is absent 
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in Matzner's (2002) treatment, which approximates with an unvarying p oc r~ 2 dependence the 
entire density profile on either side and across the He-H transition. A density drop steeper than 
r~ 2 is however expected at this core-envelope transition for radiative envelopes (Eq.|l|), and is a 
distinctive feature of the numerical pre-collapse models employed above. Using equation (^) for the 
Lorentz factor of the entire cork in an H-envelope, the additional crossing time can be much much 
smaller than the core crossing time already incurred, tu ~ r^/2cr 2 l <C t core (Meszaros h Rees, 
2001), as long as the envelope is not too extended, ru <l few x 10 12 cm, e.g. for blue supergiants 
(BSG). This is independent of whether a small fraction of the outermost cork mass is accelerated 
along the density gradient to larger Lorentz factors (Matzner & McKee 1999), as discussed in 
the next section. The only other restriction besides the envelope size is that the envelope mass 
should be smaller than that inside the stellar core, which is generally true for BSGs. Thus, the jet 
"powering" time tj needed to produce an emerging jet (and hence a successful GRB) is limited 
mainly by the core crossing time t core ~ 30 s, which is similar both for envelope-stripped CO or 
He stars and for stars as extended as blue supergiants with H-envelopes of th & few x 10 12 cm. 

The cork is composed of shocked jet and stellar plasma. The comoving width of the shocked 
jet plasma is Aj = O.20r (Meszaros &: Waxman 2001), where 9 is the asymptotic jet opening 
angle. The comoving width of the shocked stellar plasma is obtained by balancing the mass flux 
across the bow shock, 7r(9r 2 )(f3hC + v' s )(Thp) where v' s is the shock velocity in the cork frame, 
with the tangential flux of particles leaving the cylinder of (cork frame) height A s and radius 9r 
of shocked plasma, « 2i:8rA s c s p' , where c s and p' are the post-shock speed of sound and density 
respectively. Using c s = cj\[?> and p' = AT^p for the relativistic (Th 3> 1) case, which from now 
on we take as approximately valid for all cases considered, we obtain Ao — 0.29r. The Thomson 
optical depth of the shocked stellar plasma is thus 

1.5 x lO^L^/ 2 ^ 4 , H ; 
tt,s * \ 2.7 x lO^iL^V-i 4 , He ; (4) 
1.5 x l&e^Lgrtf/Z'*, CO . 



Here, 9 = 0.1#_i. The cork temperature is 

f 1-6^52 , 

T r * \l2L l ir^p 1 !l He; keV. (5) 
1 50ZV\- 1/4 P; /8 , CO. 
The corresponding (cork frame) energy carried by the cork is thus 

erg. (6) 



1 fir 1/8 r _1/ V /8 H 
i.o^ 52 r 12 5 p_ 7 , n 
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3. Jet emergence: X-ray precursor 
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3.1. Dynamics 

As the forward shock approaches the stellar edge it accelerates to relativistic velocity. 
At some finite depth d = R — r it becomes causally disconnected from the backward shock. 
This occurs at the point where during a sound crossing time of the cork, 5t ~ T^QAOr / 'c s , 
the forward shock accelerates significantly. Since oc p~ l / A (as long as the forward shock is 
causally connected to the backward one), the fractional change in the shock Lorentz factor is 
5T/T <5/9/4/9 = cStdln p/dr rj 0.4Th{c/c s )9r(n/ A)/ (R^ — r). Thus, the forward shock becomes 
disconnected at d = i?* — r ~ Tf l (d)6r, where ST /T ~ 1 (here Th(d) is the main cork Lorentz 
factor at depth d). Denoting by the value of given by Eq. (g) for r = R*/2 assuming the 
density profile of Eq. (0), we find d/R* « (r^) 4 /( n + 4 ) and T h (d) ss rt /(n+4) 0-"/( n + 4 ). Assuming 
# ~ 1/15, the value inferred by Frail et al. (2001), we find that the cork is mildly relativistic at 
the onset of "runaway" for all the progenitors (He, BSG, C/O) considered above. The dependence 
of T h (d) on 6 is not very strong, T h (d) oc p((i)" 1/4 oc (d/R)~ n / 4 oc (T n /( n+4 ) . The stellar densities 
at "runaway" onset are ~ lg/cm 3 , ~ 10~ 3 g/cm 3 and ~ 10~ 7 g/cm 3 for the C/O, He and BSG 
progenitors respectively. 

As the forward shock reaches depth d, it accelerates and breaks through the stellar edge before 
the backward shock reacts to the rapid acceleration. During jet propagation in the envelope, the 
shocked cork plasma is pushed sideways, away from the jet path, on a time scale T^Or/cs. Since 
the forward shock "runaway", from depth d outwards, takes place on a similar time scale, we 
may assume that the cork is "frozen" during this time. At break-out we will therefore have a 
cork of thickness 0A8R* with density, temperature and velocity given by the equations of §|| with 
r = R* and density p(R* — d). Ahead of this cork, there will be a faster moving shell of plasma, 
accelerated during shock "runaway" (similar to that considered in spherical supernova shocks by 
Matzner & McKee 1999). The column density of the faster plasma is smaller than that of the cork 
by a factor ~ l/(n + 1). 

During the shock runaway, the pressure behind the forward shock is smaller than 
that behind the backward shock (since the two are causally disconnected). The forward 
shock Lorentz factor is therefore smaller than that given by Eq. (||), using the density 
ahead of the forward shock: T(x) ~ (x/d) - ™/ 4 where x is the depth. The optical depth 
at depth x is t(x) ~ [tt,s/(» + ^)]{x/d) n+1 , where tt, s is the cork optical depth, implying 
T(x) ~ [tt, s /(ti + l)r(x)] n / 4( - n+1 \ Since the shocks are radiation dominated, the maximum Lorentz 
factor to which plasma is accelerated is smaller than that given by this expression for t(x) = 1. 
Using Eq. (Q), we find that the maximum Lorentz factor is of order 10 for all progenitor types 
mentioned above. 

Whether or not a fast, T ~ 10, shell is ejected ahead of the cork is unclear. As the forward 
shock accelerates and becomes disconnected from the driving jet "piston", it is likely to become 
unstable (as commonly found in the non-relativistic case, e.g. Sari, Waxman & Shvarts 2001). 
Moreover, if the progenitor plasma is magnetized, the magnetic field may also prevent large T 
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gradients. We therefore consider in the following sub-section X-ray emission from the slower, main 
cork plasma. We briefly comment here on how a fast shell ahead of the cork may affect this X-ray 
emission. 

Photons emitted by the cork would overtake the fast shell, spend some time diffusing through 
it, and then escape. Let us first consider the time delay introduced to photon arrival times seen 
by a distant observer. We show below that the fast shell does not change its radius significantly 
during the photon diffusion through it. The time delay may therefore be decomposed into two 
separate contributions: A geometric time delay due to a change in the photon propagation 
direction, which is assumed to take place instantaneously at the shell overtaking radius, and a 
diffusion time delay due to the fact that the effective velocity of the photon during the time it 
spends in the shell is smaller than c (by c/2r(rc) 2 ). 

As shown below, the time scale of photon emission from the cork is to ~ 0R*/c. Photons 
emitted radially from the cork at time t following break-out will overtake a fast shell with T(x) 
at time ti ~ 2T(x) 2 t. For t ^ to, a fast shell with T(x) > 6~ l l 2 expands significantly by t±, 
cti/R* ~ T(x) 2 e. In this case, the spread in overtaking times between photons emitted radially 
and no-radially from the cork, ~ 6 2 R*/c, is small compared to to, and we may consider geometric 
time delay of radial photons only. For an observer along the jet axis, the geometric time delay 
between two photons emitted at time t and scattered towards the observer from different points 
in the fast shell, one from a point along the line of sight and the other from the edge of the shell, 
propagating at angle 6 with respect to the line of sight (we have shown above that T(x) < 10 ^ 1/0 
so that such scattering is possible), is Q.h6 2 t\ [T(x)9] 2 t. Since T(x)6 < 1, the time spread 
introduced by this geometric time delay would not affect significantly the shape of the pulse 
emitted by the cork plasma. 

Let us now consider the diffusion time delay. For t ^ to, the optical depth of the shell 
at ti is t\(x) « r(x)(cti/ R*)~ 2 . Since the proper thickness of the fast shell is x/T(x) (due 
to the shock compression), the photon diffusion time through the shell is ~ tix/T(x)c in 
the shell frame, corresponding to St^iff ~ t\x/c in the star frame. The delay in photon 
propagation (compared to free propagation) is 5t prop « dt^if f /2T{x) 2 . For t > to we thus find 
Stprop/to ~ [T T , s /{n + l)e 2 ](t/to)- 2 T(x)- 2 ^ n+4 y n . For n = 3, the diffusion time delay is small, 
compared to to, for T(x) > [TT )S /(n + 1)# 2 ] 3 / 38 . Using Eq. (||), we conclude that this time delay 
is small for T(x) > a few. The fractional change in the shell radius during photon diffusion, 
~ dtdif f /2F(x) 2 t as 5t prop /t < 5t prop /to, is also small, justifying our assumption that the shell 
radius is approximately constant over a photon diffusion time. 

The arguments given above demonstrate that the fast shell would not affect significantly the 
pulse shape. Let us now consider its effect on the photon spectrum. The comoving post-shock 
temperature decreases towards the stellar edge, T(x) 4 oc T(x) 2 p(x) oc p(x) 1//2 , and the apparent 
temperature T(x)T oc T(x) 1 ^ 2 . The apparent temperature of the fast shell is not much larger than 
that of the cork, any expansion of the shell further reduces its temperature, and its energy content 
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is much smaller than that of the main cork. An inverse Compton contribution is also found to be 
much less energetic than other components (§ [5.2| ). Thus, we expect the fast shell to have little 
effect on the spectrum as well. 



3.2. X-ray flash at the first shock emergence 

We consider now the X-ray emission from the cork break-out. We assume that the cork 
expands at its speed of sound at the cork break-out time, c s ~ c/v3 = 0.58c in the relativistic case. 
This is the sound speed for a radiation-dominated gas, rarefaction fronts typically expanding at 
the speed of sound. As discussed before, at cork break-out the shock is at least mildly relativistic 
in the cases of interest here, so we use this limit from now on). For simplicity we replace the 
flattened cork (width = OAOr < diameter = 29r) with a sphere, and denote its radius by Ao 
= 29r/3, which has the same volume. (The spherical expansion description in the comoving 
frame is a rough but convenient approximation, which becomes increasingly justifiable as the cork 
expands beyond one e- folding of its original dimensions at the same speed in all directions.) The 
cork's radius increases (in its comoving frame) as A = Ao + c s t. We therefore define the expansion 
factor 

x = 1 + c s t/A , (7) 

so that A = Aqx, p = pqx~ 3 , T r = Tr^x -1 and tt = TTfiX~ 2 . Subscripts denote values at 
cork break-out time. Since the optical depth is very large, photons diffuse out of a thin skin 
depth at the edge of the cork. The width of the skin-depth layer out of which photons have 
escaped is tdiff — \/ctA~j3vr, where A/tt is the photon mean free path. We therefore have 
£diff/A — \/ct/3ATT = (x /TTfl) 1 ! 2 (t /to) 1 / 2 , where to = 3Ao/c. The amount of energy radiated up 
to time t is 

E rad ^^-qiaT^ = ^(^-\ 12 , (8) 



A ^/Wfl Wo 

and hence the comoving luminosity is 

L rad - J^x- 3/2 ( rV V2 • 0) 



For times t <C t exp = Aq/c s = (c/3c s )to, the expansion is not important, x ~ 1, hence L oc t~ l l 2 
and T r ~ T r ^. At later times, i > 1, L oc i -2 and T r oc t . We note that the speed of sound in 
the shocked jet plasma is always (even for (5h <C 1) relativistic, so that the rarefaction wave crosses 
the shocked jet plasma much faster than it crosses the shocked stellar plasma. Since the internal 
energy stored in the shocked jet and stellar plasma are similar (from the momentum and energy 
conservation across the shock implied by the shock jump conditions), the expansion of the shocked 
jet plasma will accelerate the expansion of the shocked stellar plasma to a velocity larger by a 
factor ~ \^2. Thus, for arbitrary (3^, whether sub-relativistic or relativistic, we can approximate 
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texp ~ to/2@h- The total energy radiated is thus 

Etot - ~^=. (10) 

One should keep in mind the potential limitations of the 1-D approximation used here, since in 
2-D numerical calculations (e.g. Aloy et al, 2000, Zhang et ah, 2002) the emergence of the reverse 
shock and rarefaction wave indicates the formation of a conical shape due to lateral expansion. 
For the parameters we used here, the 1-D approximation provides a useful description as long as 
the jet has not moved beyond much beyond ~ R*9. 

For a point on the source-observer line of sight, the characteristic break-out timescale t exp is 
shortened to t° xp = t exp /2T} l = Ao/2ThC s « 0A8r /2FhC s . The angular smearing of the light curve 
occurs on a timescale dtf° g ~ 6 2 r/2c, for which dt°^ g /t°^ p = 5T h 6c s /2c < 56T /l /2\/3. Since 9 <C 1 
and is at most a few, we have dt°^ g /t° xp <C 1 and at a distance d the flux / ~ L/d 2 of the 
observed pulse shape follows Eq. @: / oc i -1 / 2 for t < t° xp , 

2.2L 52 V_ir 12 5 p_ 7 , ri , 

t° e % =|l.4x l^L-^e^p 1 ^, He ; b, (11) 
2.4 x 10-2 L-^rfp 1 ^, CO , 

and / oc t~ 2 at later times (here d is distance of the source). The flux /q at t = t° x s p is 



approximately given by /o = r/ l (3£'o/ v / 2rr )S )/ (2t° x ! ' )ir(d/T h) in the relativistic approximation, 

f 2.2 x l^L^'l/rl^d 2 ^ H ; 
/o = 4.9 x W^L^e^/r^p^d 2 ,, He ; erg cm" 2 s" 1 . (12) 
I 2.1 x l^LTe^/rtpTd 2 ^ CO. 

The rise time of the pulse, and hence the time at which the flux peaks, is 

5.2 x KT^ 2 ^12.5, H; 
tptk « dffig « ^ 2 ^/2c = <| 1.7 x 10-^^m, He; s, (13) 

1.7 x io- 3 2 jno, CO, 



and the peak flux is 



4.4 x 10- 9 L^_ 1 /r 12 . 5 p 3 _ / 7 4 , H ; 
fpeak = { 1.4 x lO-^Ljj^-i/rnpi^, He ; erg cm" 2 s' 1 . (14) 
7.4 x W-^ 2 L% 4 9^/r w p 3 /4 , CO . 



4. Cork dynamics following break-out 



At the time t = t exp in the frame moving with Lorentz factor T^, the rarefaction wave arrives 
at the backward shock decelerating the jet; the pressure behind this shock drops, and it can no 
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longer decelerate the fast jet plasma. A new system of forward-backward shocks is formed, which 
propagates into the expanding cork plasma (see figure |l] for a schematic description). As the 
second forward shock breaks through the expanding cork, a second X-ray flash is produced. We 
consider the cork dynamics following the jet emergence in the relativistic limit, an approximation 
which is increasingly well satisfied. 

4.1. Second shock 

In the original frame of the shocked cork, moving with T^, the shocked plasma expands, 
following the penetration of the rarefaction wave, with velocity comparable to the speed of sound, 
~ c/v3. Adding (relativistically) this velocity to T^, we find that the cork expands, in the source 
frame, with 1^ « (^/S/2 + l/V2)T h = 2T h . Due to the expansion, the proper density of the 
expanding cork drops as p = pox~ 3 = AF^px -3 . Let us denote by r\ 2 the Lorentz factor (in the 
source frame) to which the rarefied cork plasma is accelerated by the second shock driven by the 
jet. As the cork density decreases with time, the second shock accelerates. Balancing the pressure 
behind the new forward shock, 3P h2 = 4{T h2 /2F' h ) 2 Ar h px- 3 c 2 = AT 2 h2 pc 2 / AT h x 3 ( assuming the 
shock is relativistic), with the pressure behind the new backward shock, 3P r = A{T j /2T] l2 ) 2 njmpC 2 , 
we find that is given by an expression similar to that determining I\, with p replaced by 
p/AT^x 3 . Using Eq. (||) we therefore find 

^ = (AT h x 3 f'\ (15) 
1 h 

Note that if we had assumed that the new forward shock driven into the cork is not relativistic, 
pressure balance gives Ph2 = x 3 ^ 2 T^ 2 /2 for the shock velocity (in the frame where cork plasma is 
at rest), which implies that the shock is indeed relativistic. 

In the frame moving with T^, the forward shock Lorentz factor is T s = y/2Tf l2 /2Tf l = (r^x 3 ) 1 / 4 . 
The shock breaks through the cork at time t for which the distance it travels, given by 
fl dt(l — l/2r^)c, equals the thickness of the expanded cork, x(i)Ao- For Th — > oo, 
T s — > oo and second shock break out occurs when xAq = c(t — t exp ) = v3Aq(x — 2), i.e. for 
x = 2v / 3/(v / 3 — 1) = 4.7. For finite T/j, crossing occurs at larger x, e.g. at x = 5 for = 4. We 
therefore adopt x = 5 for the crossing time. Note, that in this estimate we have assumed that 
is given by Eq. (^) throughout the second crossing. However, the cork is composed of shocked 
stellar plasma, and shocked jet plasma. Eq. (|l5|) holds for propagation through the shocked stellar 
plasma, while shock propagation in the shocked jet plasma is faster, due to its lower density. 
However, since the shock crossing time for = 4 (x = 5) is close to that obtained for T s — > oo, 
we neglect the small correction implied for the crossing time. 

Following the second shock crossing, the new (proper) thickness of the shocked cork plasma 
is A2 ~ xAo/[4(Tf l 2/2T' h )] = (x c /4r/j) 1//4 Ao, where x c 5 is the value of the expansion factor at 
shock break-out. The resulting second X-ray flash properties are estimated in §|5|, following the 
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same reasoning of §|3[ We give here only the estimates for the dynamical cork parameters following 
the second crossing: 

t exp2 _A 2 _ fXc y/^ 



tex P A V 4F h 

and 



, -tobs 

* " WT, = T7§k H 3 * W-^'^rli'p^, He; s. (17) 

The proper thickness A2 given by Eq. (16) is determined by the difference between the propagation 
speed of the (second) shock and the speed of the fluid behind the shock. This thickness, which is 
essentially the product of the speed difference and the shock life time, therefore grows with shock 
life time. The proper thickness saturates, i.e. ceases to increase, when the mass flux across the 
shock equals the tangential mass flux (which is proportional to the thickness). In §0, we have 
determined the proper thickness A$ of the shocked cork plasma by balancing the two mass fluxes. 
Using the same argument here would have given A2 = Ao- However, Eq. ( |l6|) gives A2 < Ao, 
which implies that the A2 does not reach its asymptotic value during shock crossing. 

Plasma expansion in the frame moving with is not isotropic. Following the first shock 
break-out, the plasma expands roughly isotropically in the frame moving with T^, to which we 
refer in what follows as the "isotropic cork frame". As the second shock propagates through 
this expanding plasma, the flow ahead of the shock is not uniform, and the geometry of the 
shock front is influenced by the shape of fluid streamlines ahead of the shock: the shock will no 
longer propagate radially in the source frame. This situation is schematically described in figure 
Since the Lorentz factor of the second shock, with respect to the fluid ahead of the shock, is 
not very large, see Eq. (|1|), we expect a large modification of the shock front geometry due to 
the non-uniform flow ahead of the shock. Assuming uniform pressure behind the shock implies 
uniform velocity of the shock front with respect to the fluid velocity ahead of the shock, i.e. radial 
shock expansion in the isotropic cork frame (moving with T/j), where the pre-shock fluid expands 
isotropically. Following the second shock crossing, therefore, the shocked cork plasma fills, in the 
isotropic cork frame, a semi-spherical shell of proper thickness A2. This shell expands radially in 
the isotropic cork frame, with a Lorentz factor Th2/2Th. The fluid propagating along the jet axis 
propagates in the source frame with Lorentz factor T/j2, and expands over an observed expansion 
time given by Eq. (|l7|). 



4 x lO-^t^-irLV- 7 ®, H; 



4.2. Third shock 

Following the second shock break-out, a second rarefaction wave propagates into the newly 
shocked cork, following which a third shock crosses the rarefied cork. We denote by the 
Lorentz factor in the source frame to which the rarefied cork plasma along the jet axis is 
accelerated by the third shock. Here too, we expect the third shock to propagate radially, 
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with Lorentz factor Tf l ^/2Tf l , in the isotropic cork frame. We define again the expansion factor 
as x = 1 + c s t /A2, where i is the proper time, measured in the fluid rest frame from second 
shock break-out. Time t measured in the isotropic cork frame is related to the proper time 
through t = (Th2/2Th)t. The radius of the expanding shell in the isotropic cork frame at 
the second shock break-out is R ~ x c c s t eX p = x c Ao ~ 26r. This radius grows with time to 
R + tc^R + {ThxlY^ct/^l = R[l + ^ffi2{T h /x c ) l lH/t exp \ = R[l + (Vs/2)t/t exp2 ] « xR. Thus, 
the proper density of the expanding cork is p = pix~ 3 = 4(r/ l2 /2r / /l )4r^/j(a;a; c ) _3 = 4Ti l 2p(xx c )^ 3 . 
Balancing again the pressure behind the new forward shock, 3P^3 = 4(r/ l3 /2r / /l2 ) 2 4r/ l 2p(^^c) _3 c 2 = 
4r 2 l3 /xi 2 / '4T h2( xx c) 3 (assuming a relativistic shock), with the pressure behind the new backward 
shock, 3P r = A(Tj/2Ti l ^) 2 njm p c 2 , we find that 1^3 is given by an expression similar to that 
determining T^, with p replaced by p/AT^ixXc) 3 . Using Eq. (|2|) we therefore find 

^ = (AT^xlx 3 ) 1 / 4 = 2 5 / 8 xi 5 / 16 x 3 / 4 rf 1 6 . (18) 

If we had assumed that the new forward shock is not relativistic, pressure balance would give 
/3/j3 = 0.5X 3 / 8 x 3 / 2 T^ 8 , which implies that the shock is relativistic. Arguments similar to those 
presented above imply that the third shock crossing also occurs at x = x c ~ 5. At the end of the 
third shock crossing we therefore have 



texp3 A3 / X c \ 



texp a V 4r^, 

and 



(19) 



t° bs 



7.5 x l^L^ 2 9- X r 2 ^p X T, H; 
Cs - = TTT^TTT; = { 7.0 x lO-SL^-ir^VJf 2 , He; b. (20) 

1.7 x w- 3 L- 2 13/32 e- ir T 6 pl 3/3 \ CO. 

Finally, we note that at the end of the third shock crossing, the radius of the expanded cork, in 
the isotropic cork frame, is x c R = x 2 Aq, and the distance 5r propagated by the cork along the jet 
axis is comparable to the stellar radius, 5r « cTj l 2X c texp2 ~ V^>x 2 T} 1 /S.q F^O-ir. 

Following the third rarefaction wave, the plasma is accelerated to T f hs ~ 2r^3 ~ 50T^^ 1G , for 
which r i /2r^ 3 pa rj/lOOr^ 1716 . Thus, following the third shock crossing the remnants of cork 
plasma in the jet path are accelerated to a high Lorentz factor. A strong fourth shock wave will 
be driven by the jet into this plasma only for Tj 3> 100. 



5. X- and 7-ray pulses from the 2nd and 3d shocks 

5.1. Thermal pulses 

The thermal emission following the second and third shock crossings can be estimated using 
arguments similar to those of §||. The energy density in the shocked cork plasma is proportional to 
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L/r 2 r| i , where subscripts i = 2,3 stand for second and third shock crossings respectively. Thus, 
T 2 = T r {T h2 /T h )~ 1 / 2 and T 3 = T^T^/Yh)- 1 / 2 , where T r is given by Eq. (§. Using Eqs. © and 
( Jj~8| ) we obtain for the temperature of the observed pulses 

6.24f>i^V 7 32 , H; 
T h2 T 2 = ^Y h ) l '^T h T r = { 2.3 x lO^fVn 716 /^ 32 , He; keV, (21) 

5.1 x lO^f 2 /^ 16 ^ 32 , CO, 

f L4xl0 1 4f 28 /r 12 3/ 5 64 p 2 V 128 , H; 
r h3 T 3 = 2 5 / 16 xf/ 32 rf 2 r h T r = 5.2 x lO^^ 128 /^ 764 ^ 1 / 128 , He; keV. (22) 

1 1.1 x l&L^/rS'"^, CO. 

Radiation emitted following the second shock crossing arrives at a distant observer lying 
along the jet axis from a cone around the jet axis of opening angle l/(Ti l2 /2Ti l ) = 2Tf l /Tf l2 
in the isotropic cork frame, where the shell radius is R = x c A$ . The angular time spread 
of the pulse is therefore given by dt°^ g2 = {2T} l /Tj l2 ) 2 R/2(2T} l )c = Y^XcAq/T^c. Comparing 
this time to the observed expansion time, t°^ s p2 = VSA 2 /2T} l2 c = V3x c ThAo/2T 2 l2 c, we find 
dt°ang 2 /t°exp 2 = 2/ a/3 1. Thus, the second pulse is "smeared" over a duration dt°^ g2 . Note that, 
had we assumed radial expansion in the source frame, the angular time spread would have been 
r/2T 2 h2 c = (r/2x c T h A )dt°^ g2 « (2.5/T h 9^x)dt°^ g2 . Thus, the resulting pulse duration is not 
sensitive to the detailed geometry of the shock front. 

The proper energy in the cork plasma following the second shock crossing is 
E 2 = 27rR 2 A 2 (aT 2 i ) = 27r(x c A ) 2 (T h x c A /T h2 )a(T h2 /T h )- 2 T r 4 = 27rA 3 (r h :r c /r /l2 ) 3 aT r 4 « 
0.3(r/ l x c /r/ l 2) 3 -E'o- Since the optical depth following the first crossing is reduced, tt, s2 = x~ 2 tt, s , 
the flux of the second pulse is f 2 « T h2 (3E 2 / ^/W^)/2tf x s p2 Tr(d/T 2 h ) » 0.3x 3 / , 

5.6 x 10-%l^/r^ 7 dli 8 , H; 
h « 0.3x 3 / = \ 6 x 10- 10 L^ /8 ^ / 1 2 /r 1 5 1 /4 p_ 3 4 8 ! He; erg cm" 2 s" 1 . (23) 
1.1 x X^L^ltlrtpJ£, CO. 

The duration of this pulse is given by Eq. (|l 

Radiation emitted following the third shock crossing arrives at a distant observer lying 
along the jet axis from a cone around the jet axis of opening angle 1/(^3/2^) = 21^/1^ 
in the isotropic cork frame, where the shell radius is R = x 2 Aq . The angular time spread 
of the pulse is therefore given by dt°^ g3 = (2T h /T h3 ) 2 R/2(2T h )c = T h x 2 c A /T 2 h3 c. Comparing 
this time to the observed expansion time, t° xp3 = v / 3A 3 /2r? i3 c = V3x 2 ThAo/2T 2 l3 c, we find 
d^angz/^expz = 2/v3 1. Thus, the third pulse is "smeared" over a duration dt°J^ g3 . If we had 
we assumed radial expansion in the source frame, the angular time spread would have been 
r/2T 2 h3 c = (r/2x 2 T h A )dt°^ g3 (2.5/x c T h 9-i)dt°^ g3 . Thus, the resulting pulse duration is again 
not sensitive to the detailed geometry of the shock front. 
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The proper energy in the cork plasma following the third shock crossing is 
E 3 = 2^R 2 A 3 (ar 3 4 ) = 2n(x 2 c A ) 2 (T h x 2 c A /T h3 )a(T h3 /T h )- 2 T r 4 = 2nAl{T h x 2 c /T h3 faT^ « 
0.3(1" \x 2 /r [ h 3 ) 3 E . Since the optical depth following the third crossing is reduced, tt, S 3 = x c 4t t,s, 
the flux of the third pulse is f 3 « T h3 (3E 3 / ^2¥^T 3 ) /2tf x s p3 Tr{d/T 2 h ) « 0.3x6/ , 

f 1.1 x Itrtlg^/r&p-rdg, H; 
/ 3 « 0.3^/ = 7.7 x lO^L^^/r^Va^ 8 , He; erg cm" 2 s" 1 . (24) 
1 1.4 x lO^I&'^/rtfpodg , CO. 

The X-ray pulses from the three shocks are schematically shown in Figure ^, appearing as 
precursors to the main burst emission. 



5.2. Non-thermal emission 

As the shock approaches the shell edge, where the optical depth ahead of the shock drops 
to ~ 1, the radiation will escape and the shock will no longer be radiation dominated. It will 
most likely become collisionless. Waxman & Loeb (2001) have shown that, for a mildly relativistic 
shock, the ratio of electro-magnetic instability growth rate and ion-ion (Coulomb) collision rate 
is i/gg/i/jj ~ 10 7 p^ . The density at the third shock crossing is smaller than the initial stellar 
density by a factor ~ x~ 6 . Thus, for all the progenitors considered here, we have vem/^u S> 1, 
and the shock is likely to become collisionless, rather than collisional, as it approaches the edge of 
the cork, where the optical depth becomes of order 1. If a fourth shock crossing takes place, i.e. 
if r,- 3> 100, the fourth shock will be collisionless as well. For the BSG progenitor, the optical 
depth of the cork shell during third shock crossing is reduced to ~ 10. In this case, therefore, a 
significant fraction of the internal energy of the cork may be emitted as radiation from electrons 
accelerated to a non-thermal distribution by the collisionless shock. 

We estimate below the non-thermal emission from the third shock as it becomes collisionless. 
In order to derive an estimate that is (nearly) independent of the details of the preceding cork 
shocks dynamics, we use the following line of arguments. The isotropic-equivalent mass of 
plasma shocked by the collisionless shock is Mql ~ 47r(r/ l 0_ir) 2 rm p / 'or, where we have taken 
rcL = T/i$-i r as the shock radius (due to expansion of the plasma during the 3rd shock crossing, 
see §||) and r ~ 1 is the optical depth ahead of the shock at the point it becomes collisionless 
{rmp/aT is the corresponding column density). We have shown that the Lorentz factors of the 
shocks, in the frame at which the plasma ahead of the shocks is at rest, are of order a few. Thus, 
the isotropic-equivalent thermal energy generated by the collisionless shock is 

f 7.7 x l^TTiT^L^e^r^pZ) 12 , H; 
E CL « TjTi M CL c 2 = I 2.6 x lO^TTiTj^^el^npZ^ 2 , He; erg, (25) 

{ 8.6 x io 44 rr i r 3 - 2 . 5 Lj/ 2 ^ irioPo x / 2 , CO, 
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where Tj ~ a few is the collisionless shock Lorentz factor (in the plasma frame). We 
have shown in §|| that the proper thickness of the shocked cork shell following each shock 
passage is close to Aq. The proper energy density behind the collisionless shock is therefore 
ucl ~ {Ecl/T j') l '^Tq L (Aot /tcl)i where tql ~ 10 is the total optical depth of the shell 
at the stage where the collisionless shock is formed. Assuming that a fraction £ e {^b) of the 
thermal energy is carried by accelerated electrons (magnetic field), the characteristic frequency 
of synchrotron photons emitted by the electrons, accelerated to a characteristic Lorentz factor 
7e ~ ie^i{m p /m e ), is 

( 7.9lf 2 e e V^UB~T h 2.^-_\ /2 r^i\ H; 
hv syn « \ 4.3 x lO 1 ^ 2 £ e 2 V^OTi, 2. 5 ezl /2 r u 1/2 , He; MeV. (26) 

[ i.4 x w 2 r 5 / 2 e e V^IJB~r j ,2. 5 ez 1 1 /2 r- 1/2 , co. 

Thus, the third thermal (hard) X-ray pulse is accompanied by a non-thermal MeV gamma-ray 
pulse, with similar energy flux and duration. 

An IC component is also expected, due to the thermal X-ray photons emitted from the cork 
which are up-scattered by the fast jet, Tj ~ 10 2 , as it emerges from the cork (Ramirez- Ruiz, et 
al. 2002). The fluence of this IC component is significantly lower than the value estimated in 
the above reference, which was based on a simplified treatment assuming that all the thermal 



X-ray photons in the cork are up-scattered by the relativistic jet. However, as shown in §3.2 



— 1/2 

only a fraction r T C 1 of the thermal X-ray photons escape the cork and are available for 
up-scattering [see Eq. (|To|)1; the value of r T ]j 2 ranges from 10~ 2 for the BSG progenitor to 10~ 4 
for the CO progenitor [see Eq. Also, as shown in §|], the fast jet emerges from the cork only 
after third shock crossing, at a time (measured in the star's frame) ^h2X c t e xp2 ~ x^Aq/c after 
the first cork break-out. This time is much longer than the transverse light crossing time of the 
jet, ~ Ao/c. Thus, most of the X-ray photons escape out of the jet's path prior to the emergence 
of the jet from the cork. We therefore conclude that the non-thermal emission is dominated by 



the collisionless emission given by Eqs. (25) and 



A second IC component may be produced if there is a T ~ 10ri fore-runner shell ejected 
ahead of the main cork, as discussed in § |3.1| , from the up-scattering of stellar photons by the 
fore-runner. The number density of stellar photons drops as r -2 , so most of the fluence is produced 
by scattering of stellar photons near R*, where i£* is the stellar radius. The (isotropic equivalent) 
energy of this IC component is ~ L^R^/cjT 2 ~ 10 37 r 2 (L*/ L & ) erg, where L* is the stellar 
luminosity. This is much smaller than other high energy components, such as equation pB]). 



6. Discussion 



We have shown that the emergence of a jet in a collapsar model of GRB leads to a series of 
successive, increasingly shorter and harder thermal X-ray pulses. These are caused by successive 
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cycles of shock waves and rarefaction waves reaching the outer and inner ends of the "cork" of 
stellar material being pushed ahead of the jet as it passes beyond the boundary of the stellar 
envelope. The Lorentz factor of the main bulk of the cork increases with each succeeding and 
increasingly relativistic shock that goes through it. A small fraction of the cork mass may be 
accelerated to larger Lorentz factors and run ahead of the cork, without changing significantly the 
development of the thermal pulses from the main bulk of the cork. 

The expansion and acceleration of the main bulk of the cork lead to a decrease of the cork 
scattering optical depth, photon diffusion times and expansion times in the observer frame for the 
successive shocks. As a result the successive shocks have increasingly higher bolometric fluences 
and shorter durations, as well as harder peak photon energies. However, for observations in a 
fixed low energy X-ray band, e.g 2 keV, the specific fluence in erg/cm 2 /keV would be reduced 
by a Ray leigh- Jeans factor (E i, s / E p k) 2 and would appear almost constant or slowly declining as 
the pulse number increases. On the other hand, in a hard X-ray detector, e.g. 20 keV, the 
successive pulses would show a steeply increasing fluence. 

The results discussed here are based on analytical 1-D calculations in the relativistic limit. 
As argued in § 2 and in § 3, the approximations inherent in such an analytical treatment are 
reasonable for the range of parameters considered, but the flows considered are quite complex, and 
definitive results must await confirmation from more accurate numerical simulations. Current 2-D 
relativistic collapsar jet simulations (e.g. Aloy et al. 2000, Zhang et al, 2002) do not so far have 
the dynamic range and resolution necessary to distinguish multiple shocks and rarefaction waves 
upon emergence of the jet, but provide valuable informations about the overall jet dynamics. The 
approximate analytical treatment presented here fulfill in the meantime an exploratory role, and 
provide insights both for future numerical and observational developments. 

The temporal structure discussed provides distinct features, which should help to identify the 
phenomenon. We have argued (§^) that since jets are likely to accelerate to significantly relativistic 
velocities as they enter a much lower density H-envelope, the crossing time for a blue supergiant 
envelope may be only slightly longer than the core crossing time t core ~ 30rn s for He or CO 
dwarf stars. The initial rise-time of the first pulse leading to the peak flux, given by the angular 
smearing time, is t° b ^ ~ [500, 20, 2] ms for a blue supergiant (H) star, a He star and a CO star, 
occurring while the cork has not yet had time to expand. This is followed by a decay / oc t~ 1 / 2 
lasting an expansion time t°^ p ~ [2, 0.15, 0.025] s for the [H, He, Co] cases. Thereafter the flux 
drops as / oc t~ 2 . The subsequent pulses have similar rise and decay times, given approximately 
by (increasingly shorter) observer-frame expansion times. 

For stars of greater compactness (smaller radius) the characteristic bolometric fluxes, timescales 
and fluences of each pulse are lower and shorter while the peak temperatures are higher. The 
typical fluences from a Hubble distance d ~ 10 28 cm range from T ~ [5 x 10~ 9 , 7 x 10~ 12 , 5 x 10~ 14 ] 
erg/cm 2 for the first pulse, to T ~ [8 x 10 _7 ,5 x 10~ 10 ,2 x 10 -12 ] erg/cm 2 for the third 
pulse, while the peak temperatures range from T°^ s ~ [1.6, 12, 50] keV for the first pulse, to 
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T pt ~ I 14 ' 52 > n °] keV for the third P ulse in I H ' He > C °] stars respectively. 

After the third shock, the cork will have reached Lorentz factors comparable to the 
observationally inferred Lorentz factors in bursts (T ~ 100). This shock may become collisionless, 
leading also to an additional non-thermal component of energy and duration comparable to the 
thermal third pulse, with characteristic photon energy £L ~ [8, 45, 140] MeV for [H, He, CO] 
stars. 

With the fluences calculated here it is apparent that it would be easier to detect the X-ray 
pulses from extended GRB progenitors, such as blue supergiant (H) stars. The more compact 
progenitors (He or CO stars) would have lower fluences and be harder to detect, except for 
nearby objects, in which case they would be distinguishable from blue supergiant progenitors 
through their shorter peak durations and harder spectra. The characteristic precursor thermal 
pulses, including the third thermal and non-thermal pulses described here, would precede the 
conventional MeV range emission by an appreciable timescale ~ t°^f p ~ 2.2L 52 1 ^ 2 O^r^ 5P-7 s ^ or 
blue supergiant (H) progenitors, or by shorter timescales ~ 0.15, 0.025 s for He or CO progenitors. 

It is conventionally assumed that the conditions for the usual successful GRB to occur is that 
the central engine powers the jet for timescales tj longer than the stellar crossing time t core ~ 30rn 
s. However, there may be cases where the jet just breaks through the outer envelope of the star 
(whether it be a He, CO or H BSG star), and does not live much longer afterwards (e.g. tj ^ 30 
s). In this case, electromagnetic detectors will see only the thermal X-ray pulses and the hard 
non-thermal gamma-ray pulse (<; 10 MeV, §[T^), which last seconds or less. These would be a new 
class of objects, perhaps related to X-ray flash bursts and/or to short GRBs, especially if such 
events are more common, and their distance is smaller. 

The characteristic timescales as well as the ratios of the amplitudes and hardnesses of the 
pulses differ substantially, depending on the density and extent of the stellar envelope, which 
provides a potentially valuable diagnostic for pinning down the conditions in the pre-burst stellar 
progenitor. Instruments to be flown in the next several years, such as Swift, GLAST, AGILE and 
others, should be able to detect such precursor signals in some bursts, and provide valuable clues 
concerning the GRB mechanism and their progenitors. 

This research has been supported by NSF AST0098416, BSF 9800343, Universities Planning 
& Budget Committee, Israel, and NASA NAG5-9192. 
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Fig. 1. — The top panel presents a schematic description of the double shock structure at the head 
of the jet, as seen in the frame where the shocked plasma is at rest, during jet propagation within the 
star. The dot-dashed line shows the pressure distribution. The bottom panel presents a schematic 
description of the flow after forward-shock break-out. A rarefaction wave propagates backward 
into the shocked plasma, and the plasma ahead of this wave expands with velocity comparable to 
the speed of sound. When the rarefaction wave reaches the backward shock, a new double-shock 
structure similar to that shown on the top panel is formed, propagating through the expanding 
shocked plasma. 
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Fig. 2. — A schematic description of the flow during second shock propagation through the cork, 
as seen in the "isotropic cork frame" moving with T^. 
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Fig. 3. — X-ray precursors of increasing flux and hardness and decreasing duration, resulting from 
successive shocks and rarefaction fronts moving across the cork of stellar material ejected by the 
GRB jet. 



